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Abstract: In this work we consider the SU(N) gauge theory with two Dirac fermions 
in the adjoint representation, in the limit of large N. In this limit the infinite-volume 
physics of this model can be studied by means of the corresponding twisted reduced model 
defined on a single site lattice. Making use of this strategy we study the reduced model 
for various values of N up to 289. By analyzing the eigenvalue distribution of the adjoint 
Dirac operator we test the conformality of the theory and extract the corresponding mass 
anomalous dimension. 
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1. Introduction 

The SU{2) gauge theory with two adjoint Dirac fermions, known as Minimal Walking 
Technicolor (MWT) [1, 2], has been the subject of many lattice studies, all of which have 
found it to be a conformal theory with a fairly small mass anomalous dimension 7 * [3, 4, 
5, 6 , 7]. The most recent measurement obtained by fitting the mode number of the Dirac 
operator gave a very precise value [ 8 ]. The mode number method has also been used to 
follow the running of 7 over a range of energy scales for the SU (3) theory with many light 
fundamental fermions [9]. 

The large N version of MWT, the SU{N) gauge theory with two adjoint fermions, is 
interesting for several reasons. From a phenomenological point of view, it is expected to be 
similar to the SU{2) theory. For example, the universal Hrst two perturbative coefficients 
of the beta function are independent of N. Hence, as in the case of MWT, they point 
towards the existence of an infrared fixed point with a mass anomalous dimension that 
is also independent of N. Moreover, numerical results for the mass anomalous dimension 
for the SU(2) and SU(3) theories appear to be in agreement [10], suggesting that this N- 
independence may be a good approximation all the way down to = 2. From a more 
theoretical point of view, the large N theory is better suited for connecting with results 
obtained from different approaches, such as the AdS/CFT correspondence. Fortunately, 
the numerical study of the infinite volume theory at large N is made possible by the 
concept of large N volume independence. This implies the equivalence with a single site 
lattice reduced model, for which simulations can be performed at large values of N, that 
would be prohibitively expensive on a conventional lattice. In this context, the study of 
large N Yang-Mills theory with adjoint fermions has attracted much attention [11]- [20]. 
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In this work we will be using the twisted reduction technique [21, 22], For the adjoint 
fermion case the specific form of the action has been given in Ref. [23]. This model has 
been shown to lead to a softer N dependence than the Adjoint Eguchi-Kawai model with 
periodic boundary conditions [12, 23]. The twisted model depends on the choice of the 
twist tensor. Here we will follow the same symmetric twist prescription as for the pure 
gauge theory, in which N is taken as the square of an integer number N = L?', and the 
flux through each plane is equal to ±.kL (modulo N). With appropriate values of the 
integer k, this choice has proven effective in avoiding symmetry breaking for the pure 
gauge theory [24]. An important advantage of the twisted reduction method is that the 
dominant 1/A^ corrections amount to finite size effects on an lattice. This allows an 
estimate of the values of N at which the simulations should be performed. In this work we 
will be using values of N up to 289, corresponding to lattices of size 17^. 

In summary, the purpose of this paper is to analyze the behaviour of the SU{N) 
gauge theory with two flavours of adjoint fermions in the large N limit. Our main goal 
is to determine whether the theory has indeed a non-trivial infrared fixed point (IRFP) 
and to measure the mass anomalous dimension at this fixed point. In previous papers 
some of the present authors studied the behaviour of Wilson loops and the corresponding 
string tension [25, 26, 27, 28]. Although the results were consistent with the conformal 
behaviour characteristic of an IRFP, the extraction of the mass anomalous dimension had 
large uncertainties. Our methodology here will be based on an alternative procedure which 
has produced very precise estimates in the study of MWT [ 8 ]. Preliminary results have 
been presented in Refs. [29], [30]. 

The strategy is to determine the anomalous dimension from the structure of the eigen¬ 
value density p{ijj) of the massless Dirac operator Ip. The eigenvalue density is defined 
as 


p{uj) 


lim — 5{w — 

V^oo V 


Wk) 


(I.l) 


k 


where the sum runs over all eigenvalues iojk oi Ip. In a mass-deformed conformal field 
theory (mCFT), this quantity should vanish for u) —>• 0 as [31] 


lim lim p{uj) (x , ( 1 - 2 ) 

m—^O V^oo 


where 7 * is the mass anomalous dimension at the infrared hxed point, V is the lattice 
volume, and m is the mass. This behaviour should be contrasted with the one characteristic 
of a chirally broken theory where the eigenvalue density does not vanish at the origin. 

In this paper we will use the previous idea to determine 7 * for the large N gauge 
theory with two adjoint quarks. It is clear from the previous formula that it is crucial to 
work in the region of very small masses and keeping finite volume effects under control. 
As a reference we will compare our result with those obtained for the pure gauge theory 
{uf = 0 ), for which the eigenvalue density has the characteristic behaviour of a chirally 
broken gauge theory. 

The structure of the paper is as follows. In the next section we will collect all the 
technical aspects concerning the simulation and the extraction of the mass anomalous 
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dimension from the data. In the following we will present the results of our analysis. The 
paper ends with the presentation of our conclusions. 

2. Methodology 

As explained in the introduction, our approach to the large N limit is based on reduction. 
Hence, we simulate the twisted reduced SU(N) model on a single site with two adjoint Dirac 
fermions [23]. In the large N limit the theory is equivalent to the infinite volume lattice 
gauge theory. For finite N = the corrections amount to finite volume corrections in an 
lattice. Thus, it is important to keep track of the A^-dependence which translates into 
the equivalent finite volume corrections. For that purpose we have performed simulations 
at values of N ranging from 16 up to 289, the latter corresponding to an effective lattice 
volume of 17^. Our study has been done at two values of 6 , 0.35 and 0.36, and a large 
number of n values. The number of configurations used for the calculation of the eigenvalue 
spectrum at each value of b, k and N are listed in Tab. 1. In addition we calculated the 
eigenvalue spectrum of the Uf = 0 theory at 6 = 0.35, 0.36. For N = 121, 289, we used 10 
conhgurations with k = 0.170,0.175, 0.180, 0.185,0.190, and for N = 841, b = 0.36 we used 
4 conhgurations for k = 0.190. 




N = 16 

II 

to 

II 

N = 121 

N = 289 


k 

1 

2 

3 

3 

5 


k/VN 

0.25 

0.40 

0.29 

0.36 

0.41 

b 

K 

N = 16 

II 

to 

II 

N = 121 

N = 289 

0.36 

0.160 

20k 

1200 

500 

20 

20 


0.165 

- 

- 

- 

20 

19 


0.170 

- 

- 

- 

20 

15 

0.35 

0.160 

- 

- 

- 

40 

20 


0.165 

- 

- 

- 

20 

20 


0.170 

- 

- 

- 

20 

13 


Table 1: Number of conhgurations used for calculating the eigenvalue spectrum at each 
value of b, k, N and hux k. The integer parameter k satishes kk = 1 (mod ^/N). All 
conhgurations are separated by at least 25 molecular dynamics updates, and at N = 289 
they are separated by 125 molecular dynamics updates. 

As explained in the introduction, we have chosen the symmetric twist conhguration, 
with values of the hux integer parameter k given in Tab. 1. These fulhll the condition 
k/y/N >1/9 which was found necessary for the pure gauge theory (TEK model) to respect 
the center symmetry [24] . This symmetry is a necessary ingredient in the proof of reduction 
by Eguchi and Kawai. The addition of light adjoint fermions should help in preserving the 
symmetry but, as shown in Ref. [23], adhering to the condition allows the study of the full 
range of k values and leads to a smoother N dependence. As an example of the behaviour 
of Polyakov loops, which act as order parameters of the center symmetry, in Fig. 1 we 
display the expectation value of the modulus of the unit winding loop as a function of N. 
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By definition, this quantity is always positive, but as seen in the figure its size decreases 
with N for all values of n. 


Polyakov Loop vs 1/N, nf=2, b=0.36. 



limit for reduction to hold, which it does for all values of k, where k = 0 corresponds to 
quenched Uf = 0 data. 

For each configuration we compute the low-lying spectrum of the modulus square of 
the massive lattice Wilson Dirac operator in the adjoint representation. This operator is 
positive definite and in the naive continuum limit corresponds to a?{—lj) +m^). Thus, its 
eigenvalues, labelled (oD)^, are related to those of ij) by the expression 

auj = \/ (aD)2 — rrfla^ , (2.1) 

where a is the lattice spacing and m is the adjoint quark mass. The lowest eigenvalue of 
our lattice operator defines the spectral gap. In the continuum it is bounded from below 
by a?rin?. In the case of QCD and for the lattice Wilson Dirac operator, the median of the 
gap distribution was found empirically [32] to satisfy the relation 

aDo oc (2.2) 

so that oDo is proportional to am. In our case, however, we expect that the bound is not 
saturated at finite N. The reason being the absence of zero-momentum quark states in the 
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reduced model. Hence, quarks are then produced with at least the minimum momentum 
27: jL^ whose square decreases linearly with 1/N. 

Since our main goal is the determination of the mass anomalous dimension with the 
idea presented in the introduction, it is important to study the region close to the critical 
point and for large values of N to minimize small effective volume corrections. For that 
reason, our main analysis was based on the study of the lowest 2000 eigenvalues (aH)^ 
at = 289 and the lowest 1000 eigenvalues at = 121. The set of values of b and k 
were given in Tab. 1. The computational cost increases considerably as we approach Kc, 
explaining the smaller number of configurations for that case. Fortunately, the distribution 
of the lowest lying spectrum does not seem to fluctuate strongly at those values of N. 

In determining the value of 7 * from the distribution of eigenvalues there are certain 
alternative procedures which we will describe below. 


2.1 Determination of the mass anomalous dimension 7 * 

2.1.1 Determining 7 * from a fit to the spectral density 

In the continuum 7 * could be determined by fitting the spectral density to the form expected 
for a mass-deformed conformal field theory, Eq. (1.2). However, in order to compare to the 
lattice data, it is more convenient to look at the spectral density of —Ip , given by: 

= lim - wl) oc (u;2)t+^ . (2.3) 

V^oc V ^^ 
k 


This quantity is obtained on the lattice by counting the number of eigenvalues of the modu¬ 
lus square of the Wilson Dirac operator within a bin of size around Representing 

this number by AA(aH, A), the lattice spectral density is given by: 

PL((anf) = j^Ar(an,A), (2.4) 

where represents the lattice volume on the reduced lattice. The continuum formula for 
p(cij^) gives the following parameterisation for the lattice data: 

Pl((oH)^) = B [(aH)^ — (am)^] , (2.5) 


allowing the determination of the three free parameters: B, {am)'^ and 7 *. The lowest 
part of the eigenvalue distribution is the one most affected by finite volume and finite 
mass effects, hence the fits have to be performed in an intermediate range of eigenvalues 
a^min < o,VL < aQmax, which preserves the separation of scales on the lattice, 


<C am oH <C 1 . 

Vn 


( 2 . 6 ) 


From a practical viewpoint the 3-parameter fit demands very precise data in a wide range 
of eigenvalues and induces strong correlations between the parameters. In some cases 
the number of parameters can be reduced by assuming that the mass is negligibly small 
(am = 0 ). Alternatively we can continue to work at finite mass but use information coming 
from the smallest eigenvalue to fix the parameter (am)^ of the fit. 
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(a^2o)^ vs 1/N for b=0.36. Fit to (am)^ + c/N 



Figure 2: Median of the spectral gap of the modulus square of the Wilson Dirac operator, 
(aDo)^, vs 1/N. The htting form (aDo)^ = (am)'^ + c/N fits all the data well. 

2.1.2 Determining 7 * from a fit to the mode number 

Alternatively we can follow the procedure introduced in Ref. [ 8 ] and extract 7 * from the 
mode number i^(D) of the Dirac operator. It is simply dehned as the number of eigenvalues 
of {—l/> + m^) below some value Hence, it is given by V times the integral of the 
eigenvalue density. We can split this integral into two parts as follows 

z^(D) = 2V / p{uj) dio + 2V / p{u}) duj, (2.7) 

Jo j yj 

The first part contains the range of eigenvalues which are more sensitive to finite volume 
and/or finite mass effects. For the second part we can insert Eq. (1.2) and perform the 
integration to give 

u{Q) ~ — A [lljji — ^+T'* + A — m^] ^+t'* . (2 .8) 

To determine i/(D) on the lattice we just simply count the number of eigenvalues of the 
modulus square of the lattice Wilson massive Dirac operator below some value (aD)^. The 
continuum formula for the mode number implies that the lattice data can be parametrized 
as follows: ^ 

UL{aVt) ~ 1^0 + ^ [(oD)^ — (am)^] , (2.9) 

where 

^0 = — A [(oD//?)^ — (am)^] ^+t'* . (2-10) 


- 6 - 
















Eq. (2.9), normalised dividing by the lattice volume, is the expression used in Ref. [8] to 
ht the lattice data. This allows the determination of its four free parameters (z/q. A, 
and 7 *). 

For the same reason as discussed before for the spectral density, one can attempt to 
reduce the number of parameters in the fit. If we assume that finite volume (i.e. hnite 
N) effects are negligible we might take close to threshold making uq negligibly small 
compared to the other term. This leads to a simplihed expression 

z^L(an) ~ A — (am)^] . (2T1) 


which can be used to fit its three parameters {A,am,'yAj to the modenumber data in a 
range of eigenvalues G satisfying; 


-^= <C am <C aQ.iR < aft < ailuv ^ 1- 
VN 


( 2 . 12 ) 


If we set am = 0 in Eq. (2.11) we can reduce the number of free parameters even 
further. This is for example the strategy adopted in Ref. [9] where a 2-parameter fit to 
the mode number is used, 

~ ^[(all)^] i+T"*. (2-13) 

The sensitivity to the ht function, the volume, the mass parameter, or the htting range 
will be used to estimate the systematic error in the determination of the mass anomalous 
dimension. 


b 

K 

N = 121 

N = 289 

N = oo 

0.36 

0.160 

0.0803(6) 

0.0585(3) 

0.0429(7) 


0.165 

0.0433(4) 

0.0224(2) 

0.0074(5) 


0.170 

0.0279(4) 

0.0096(2) 

-0.0036(4) 

0.35 

0.160 

0.0997(9) 

0.0815(4) 

0.0683(9) 


0.165 

0.0530(5) 

0.0346(2) 

0.0214(6) 


0.170 

0.0281(4) 

0.0105(4) 

-0.0021(7) 


Table 2: Lowest eigenvalues squared for each b, k and N, along with the extrapolation to 
N = oo which corresponds to the mass parameter (am)^. 


3. Results 

3.1 Analysis of the two flavor case 

In this section we will present the results of our study. Our hrst step is the analysis of the 
spectral gap of the hermitian Wilson-Dirac operator. To study how this quantity behaves 
as a function of N and zt, we measure the lowest eigenvalue at b = 0.36 on a range of 
conhgurations for N = 25 — 289 and k = 0.130 — 0.170. As argued in the previous section, 
in the twisted model we expect the median of the gap distribution (aflo)^ to differ from 
(am)^ by a finite volume correction which, interpreted as non-zero momentum contribution. 
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b=0.35 


b=0.36 




Figure 3: Probability distribution of the spectral gap of the Wilson Dirac operator for 
N = 289 and various values of k. The lines are fits to a gaussian distribution, Eq. (3.1). 



K 


Figure 4: We display N^a for our ensemble of nj = 2 configurations for b = 0.36 and 
b = 0.35 (displaced by 3 vertically for clarity), where is the width of the probability 
distribution of the spectral gap of the Wilson Dirac operator. 


should depend linearly on 1/N. Indeed, a linear fit of this kind seems to describe our non- 
perturbative data quite well, as shown in Fig. 2. This allows us to determine the mass 
parameter (am)^ in the large N limit, listed in Tab. 2. 

We have analysed the probability distribution of the spectral gap p(aDo) of the Wilson 
Dirac operator. In the large volume limit of QCD, the analysis in Ref. [32] showed that 
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Histogram of eigenvalue density: b=0.36,K=0.16 



Figure 5: Spectral density pL{{aQ)‘^) for b = 0.36, k = 0.160, with a bin-size = 0.01. 
There is agreement at all eigenvalues (except the very lowest couple of bins) between N = 
121 and N = 289. For N = 49 there is a qualitative difference for the lower eigenvalues, 
but the large ones are also in agreement. 


the distribution is gaussian: 

p(ano) oc exp | — 


2(e«) 


(afio - I 


(3.1) 


with median proportional to the bare current quark mass and width scaling with the volume 
and the lattice spacing approximately as = o? jyfV. In the reduced lattice this would 
imply = 1/A. Fig. 3 shows p{aQo) for N = 289 at 6 = 0.35 and 0.36 and several values 
of K. A gaussian fit describes the data well. The htted distribution widths multiplied by 
N are displayed in Fig. 4 for our Uf = 2 configurations at several values of N and k. Our 
results follow rather well the behaviour also reported in QCD. 

Let us now move on to describe our results for the distribution of eigenvalues. One 
of the main points is to analyze the A-dependence of this distribution. We already saw 
that this dependence affects the gap of the spectrum, but we expect this effect to have a 
small impact for higher eigenvalues. This can be seen in Fig. 5, which shows a histogram 
of the number of eigenvalues as a function of (afl)^, for b = 0.36, k = 0.160. Comparing 
different values of A, we see agreement between A = 289 and A = 121 in all but the 
first two bins. For A = 49, for (all)^ < 0.4 the behaviour is qualitatively different, but 
for higher eigenvalues we again see agreement with larger values of A. For A = 25 and 
A = 16 there are strong oscillations in the distributions, which are presumably the sum of 
the distributions of individual eigenvalues with the allowed discrete momenta. 
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Histogram of eigenvalue density: nf=2, b=0.35, 7=0.27 Histogram of eigenvalue density: np2, b=0.36, 7=0.27 




Figure 6 : Eigenvalue density distribution at 6 = 0.35 (left) and b = 0.36 (right) for k = 
0.165,0.17 and N = 121,289. 


We have seen that, at least for the two biggest values of N, the eigenvalue distribution 
roughly coincides beyond a certain threshold value. The question now is to see if this 
distribution behaves as expected from the IRFP hypothesis and to extract 7 * from it. In 
the previous section we gave two alternative methods of fitting the data. One is to compare 
the eigenvalue distribution with Eq. (2.5). The other is to compute the mode number and 
fit it to Eq. (2.8) or to its simplified expression Eq. (2.11). 

In performing a fit one has to select the range of values fimin < < Umax to be fitted. 

A lower value of fimin increases the sensitivity to the value of the mass (om)^ but also risks 
to be more affected by finite effective volume (hnite N) corrections. For the mode number 
fit the same is expected to happen for the parameter vq. Furthermore, the narrower the 
fitting range the stronger the correlations among parameters leading to high uncertainties 
in 7 *. 

If we use only the data which is least affected by finite volume and finite mass effects, 
we can produce our most precise determination of 7 *: Hence, we take our data of = 289 
and K = 0.17 and ht the distribution of eigenvalues to Eq. (2.5) with set to zero. The 

upper edge of the htting range (allmax)^ covers almost all of the N = 289 data. For the 
lower edge (oHmin)^ the cut is set to twice the lowest eigenvalue for N = 121. The result 
for b = 0.36 is 7 * = 0.268(2) and for b = 0.35 is 7 * = 0.271(1). It is remarkable that both 
values of b give consistent results within the purely statistical 1% errors. To give an idea of 
the quality of the fit we display the data in Fig. 6 together with the best fit. In the figure 
we also include the data of = 121 at both values of b. The fitted function also describes 
well the behaviour of the data at A^ = 121 , except for the smallest eigenvalues where finite 
volume effects should be mostly felt. The continuous line going through the k = 0.165 data 
was obtained fitting the corresponding data with 7 * fixed to the value obtained at k = 0.17 
and the mass square (am)^ to the large N extrapolated value given in Tab. 2. 

The analysis of the previous figure and fits indicates that our data look consistent with 
the predictions of an infrared fixed point with mass anomalous dimension close to 0.27. 
In order to substantiate the claim of conformality, it is important to check the systematic 
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Fitted Y vs lower bound of fit range, nj=2, k=0.170. 



Figure 7: Fitted 7 * as a function of (aOmin)^ from fitting the spectral density data to 
Eq. (2.5). For each point, a range of masses (om)^ between zero and the lowest eigenvalue 
at = 289 are used, and the upper fit range bound (allmax)^ is also varied over a range 
of values from 0.25 up to the largest value for which we have data. The black line shows 
our best estimate for 7 * with statistical error only. 


uncertainties involved in our analysis and to compare our results with those corresponding 
to theories which are not conformal in the IR. The latter will be done in the next section 
where we will compare our n/ = 2 results with those obtained for n/ = 0. What we 
will now present is an evaluation of the systematic errors. This point will be analysed by 
estimating the effect of finite N corrections, finite mass corrections and sensitivity to the 
fitting range. We will also test what results come out if we use the mode number instead 
of the eigenvalue density. 

A good summary of the effect of all systematics on the value of 7 * is provided by 
Fig. 7. Here we display different determinations of 7 * for the N = 289 and N = 121 data 
by varying Umax, and The latter is varied within the range extending from 

zero to the minimum N = 289 eigenvalue. The data is displayed as a function of (oHmin)^- 
For a large range of x-axis values all determinations of 7 * (obtained by using different 
values of Hmax and (am)^) fall within a horizontal strip whose width serves as an upper 
bound to our systematic error 57 * = 0.05. For larger values of (oHmin)^ the fitting range 
narrows, obviously leading to a wider spread of values of 7 *. 

To better analyze the dependence of the fitted value of 7 * with the remaining param¬ 
eters, we display in Fig. 8 its correlation with the mass parameter used in the fits for all 
values (oHmax)^ and (aflmin)^ < 0.07. This shows that the main source of systematic errors 
is indeed the value of the mass. 

A similar analysis can be done for the mode number. Results are essentially compatible. 
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Fitted y and (am)^. n,=2, N=289, k=0.170. 



b=0.36 I — I- 
b=0.35 


03 -0.02 


- 0.01 


0.01 


0.02 


0.03 


0.04 


(amf 


Figure 8 : Fitted values of {am)"^ and 7 * when both are taken as free parameters when 
fitting the spectral density data to Eq. (2.5), using N = 289 configurations. There is a 
clear correlation between the two parameters. 


Mode number of eigenvalue density: r\f=2, b=0.35, y=0.27 



Mode number of eigenvalue density: 0^=2, b=0.36,7=0.27 



Figure 9: Fits to the mode number atb = 0.35 (left) and b = 0.36 (right) for k = 0.165, 0.17 
and N = 121,289. 


In this case, one has an issue about which fitting formula one should use. Although the 
4-parameter formula Eq. (2.8) should apply in all regions, there are strong correlations 
between the fitted parameters, so that for example a large range of values of 7 * can produce 
a good fit by a suitable choice of the other 3 ht parameters. The more restrictive fitting 
formulas better constrain the htted quantities, but have a limited range of applicability. For 
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N=289, b=0.36 


Nf=2 


N,=0 




Figure 10: Eigenvalue density distribution, in log-log scale, at 6 = 0.36 for nj = 2 and 
re/ = 0. The lines correspond to the fits to Eq. (2.5) described in the text. 

example Fig. 9 shows a fit to Eq. (2.11) using the same configurations and fit ranges as for 
the eigenvalue density in Fig. 6 , which gives 7 * = 0.270(2) for b = 0.36 and 7 * = 0.272(1) 
for b = 0.35, which are in good agreement with the numbers determined from the eigenvalue 
density. 

3.2 Comparison between zero and two flavors 

The comparison between our re/ = 2 results with those obtained for re/ = 0 is essential 
to substantiate the claim of conformality. The latter theory is not conformal and should 
display a different behaviour. Indeed, for the pure Yang-Mills theory it makes no sense to 
speak of 7 * itself, since there is no infrared fixed point. However, it still makes sense to 
study the behaviour of the spectral density and mode number distribution as a function of 
its argument. 

To carry on the previous study, we generated pure gauge configurations at the same 
values of b and N and a range of k values, and repeated our previous analysis for this data. 
The K values were chosen in such a way as to explore similar small values of the minimum 
eigenvalue. 

The fits to the spectral density for the largest values of re were qualitatively as good 
as those of re/ = 2. Fig. 10 shows, in log-log scale, the fits to the lattice spectral densities 
for re/ = 2 and re/ = 0 at 6 = 0.36. The values of 7 * in this plot are obtained by fitting the 
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Figure 11: a) Dependence of the extracted value of 7 * with the mass (am)^ for nf= Oand 
nf = 2. b) Dependence of the extracted value of 7 * with {m)‘^/a for nj= Q 


spectral densities to Eq. (2.5). The mass (am)'^ is varied in the range 0.5(aDo)^ < < 

(aDo)^, the number of bins is varied from 40 to 80, and the lower edge of the fit range is 
varied in the range 1.5(aDo)^ < < 2(aDo)^. The final point is an average over all of 

these choices for many bootstrap replicas of the data. Keeping 7 * fixed to this average, the 
fits are repeated for each set of data in order to determine am. For this final fit the lower 
edge of the ht range is set to 2(aDo)^. In all cases the degree of freedom of the fit is 

below 2. Hence, from the quality of the ht to the spectral density one cannot deduce the 
presence of an infrared hxed point. 

One remarkable difference appears when looking at the dependence of the extracted 
value of 7 * with the mass (am)^. The result is shown in Fig. 11a. The value of 7 * seems 
quite stable for the nf = 2 case. This is what one expects in the vicinity of an infrared 
hxed point as the anomalous dimensions tends to a constant at the hxed point. The result 
for Uf = 0 is quite different, showing a pronounced drop as we move away from the critical 
value of K. For the smallest masses the htted value of 7 * reaches as high values as 0.8. 
Notice that a value equal to 3 would imply a constant value of the spectral density at 
the origin. Our data show a growing 7 * for lighter masses but do not yet reach the value 
of 3 predicted by chiral symmetry breaking, through the Banks-Casher formula. This is 
probably due to hnite volume and/or hnite mass effects. 
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Another marked difference between the two cases appears in the dependence of the 
spectral densities on the bare inverse coupling b. In the conformal case one expects the 
coupling to be a marginally irrelevant operator close to the IRFP, in contrast to the nj = 0 
case where it is marginally relevant and determines the lattice spacing. The values of 7 * 
displayed in Fig. 11a for nj = 2 do indeed have a very small dependence on b. However, this 
dependence is large for the case of nj = 0. In the quenched case, the results corresponding 
to the two different couplings only show scaling when expressed in physical units in terms 
of the string tension. This is shown in Fig. 11b, where we have used the data for the 
re/ = 0 string tension obtained in Refs. [33] {aa^ = 0.09 and crre^ = 0.043 for b = 0.35 and 
b = 0.36 respectively). 

4. Conclusions 

We have performed a measurement of the mass anomalous dimension 7 * of the SU(A^) 
gauge theory with two adjoint Dirac fermions, in the large N limit using the concept of 
large N twisted reduction. Results from a single site lattice model at large values of N, 
have the expected qualitative behaviour of the spectral density and mode number of the 
adjoint massless Dirac operator. The distribution for small masses (extracted from the 
lowest eigenvalue) can be well-fitted with the expectations of an IRFP. From the data we 
extract a value of 7 * = 0.269±0.002±0.05, where the first error is statistical and the second 
one systematic. This value is similar to previous lattice determinations of this quantity for 
the SU(2) theory. 

Does our result provide conclusive evidence of the presence of an infrared hxed point 
for the SU(oo) gauge theory with 2 flavours of adjoint fermions? To try to answer this 
question we repeated the analysis for the re/ = 0 case, which is known to have a completely 
different behaviour at criticality. However, we observe that the spectral densities at fixed 
b and small quark mass, can also be fitted with the same formulas with a larger value of 
7 *. This conveys a word of warning about drawing conclusions about the existence of an 
IRFP only from the capacity to fit the spectral density or mode number with a powerlike 
distribution. Nevertheless, there are marked differences between the behaviour observed 
in the re/ = 0 and re/ = 2 cases. One of them is the dependence of the result on the bare 
coupling b. The re/ = 2 results for our two values of b, 0.35 and 0.36, are consistent with 
each other. This is not the case for the re/ = 0 data. Although insensitivity to the bare 
coupling b is certainly the expected result for an IRFP, it is difficult to exclude the fact that 
this is not simply due to the smaller value of the beta function when adding fermions in the 
adjoint. The second difference refers to the change of behaviour as we approach criticality 
K —>• Kc- The extracted value of 7 * for re/ = 2 remains fairly stable as one expects if the 
behaviour is indeed dictated by the presence of an IRFP. On the contrary for the re/ = 0 
data we observe a pronounced rise of the value of 7 * as we approach criticality. 

To improve on these results using this method, smaller fermion masses and larger 
volumes (i.e. larger values of N) would be required. Since finite volume effects are the 
dominant limitation, one possibility for future work would be to extend the single site 
lattice to a 2 ^ lattice. 
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